Four-body Efimov effect 
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We study three same spin state fermions of mass M interacting with a distinguishable particle of 
mass m in the unitary limit where the interaction has a zero range and an infinite s-wave scattering 
length. We predict an interval of mass ratio 13.384 < M/m < 13.607 where there exists a purely four- 
body Efimov effect, leading to the occurrence of weakly bound tetramers without Efimov trimers. 
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In a system of interacting particles, the unitary limit 
corresponds to a zero range s-wave interaction with infi- 
nite scattering length [l|. In particular, this excludes any 
finite energy two-body bound state. Interestingly, in the 
three-body problem, the Efimov efi'ect may take place , 
leading to the occurrence of an infinite number of three- 
body bound states, with an accumulation point in the 
spectrum at zero energy. This effect occurs in a variety 
of situations, the historical one being the case of three 
bosons, as recently studied in a series of remarkable ex- 
periments with cold atoms close to a Feshbach resonance 
[3| . It can also occur in a system of two same spin state 
fermions of mass M and a particle of another species of 
mass m, in which case the fermions only interact with the 
third particle, with an infinite s-wave scattering length: 
An infinite number of arbitrarily weakly bound trimers 
then appears in this 2 + 1 fermionic problem if the mass 
ratio a = M/m is larger than ac{2; 1) ~ 13.607 0]. 

The four-body problem has recently attracted a lot of 
interest [1] . The question of the existence of a four-body 
Efimov effect is however to our knowledge still open. We 
give a positive answer to this question, by investigating 
the 3 -|- 1 fermionic problem in the unitary limit. We 
explicitly solve Schrodinger's equation in the zero range 
model [4I and we determine the critical mass ratio to have 
a purely four-body Efimov effect in this system, that is 
without Efimov trimers. 

In the zero-range model, the Hamiltonian reduces to a 
non-interacting form, here in free space 
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with mi = TO2 = 'TIS = M and 7714 = m. The interac- 
tions are indeed replaced by contact conditions on the 
wavefunction, ■(/'(ri, r2, ra, r4), where r^, i = 1, 2, 3 is the 
position of a fermion and r4 is the position of the other 
species particle: At the unitary limit, for i — 1,2,3, there 
exist functions A, such that 
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when Ti tends to for a fixed value of the i-4 centroid 
Ri4 = {Mri + mr4)/{m + M) different from the positions 
of the remaining particles r/j, /c 7^ i, 4. The wavefunction 
is also subject to the fermionic exchange symmetry with 
respect to the first three variables r^, i = 1,2,3. 

In what follows, we shall assume that there is no 
three-body Efimov effect, a condition that is satisfied 
by imposing M/m < q;c(2; 1) — 13.607. The eigen- 
value problem Hip = Et/j with the contact conditions 
in Eq.© is then separable in hyperspherical coordi- 
nates [5|. After having separated out the center of 
mass C of the system, one introduces the hyperradius 

r . -,1/2 

R = J2t=i - C)^/m , with m = (3A/ + to)/4 

the average mass, and a set of here 8 hyperangles whose 
expression is not required. For a center of mass at rest, 
the wavefunction may be taken of the form 



V(ri,r2,r3,r4) = 
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/(ri) is given by the solution of a Laplacian eigenvalue 
problem on the unit sphere of dimension 8, which is non 
trivial because of the contact conditions. On the con- 
trary, the hyperradial part F is not directly affected by 
the contact conditions, due iiiparticular to their invari- 
ance by the scaling Xri l6|, and solves the effective 
2D Schrodinger equation 
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The quantity is given by the hyperangular eigenvalue 
problem. It belongs to a infinite discrete set and is real 
since there is no Efimov effect on the unit sphere (i? ^ 0), 
that is here no three-body Efimov effect. 

Mathematically, Eq.Q admits for all energies E two 
linearly independent solutions, respectively behaving as 
R^' for i? 0. If s^ > 0, one imposes F{R) - R" , 
with s > 0, which is correct except for accidental, non- 
universal four-body resonances (see note [43] in (H), and 
Eq.Q then does not support any bound state. On the 
contrary, if s^ < 0, in which case we set s = iS. S > 0, F 
experiences an effective four-body attraction, with a fall 
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to the center leading to a unphysical continuous spectrum 
of bound states To make the model self-adjoint, one 
then imposes an extra contact condition [3], as in the 
usual three-body Efimov case Q : 
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where the four-body parameter Rf depends on the mi- 
croscopic details of the true, finite range interaction @. 
With the extra condition Eq.(IS]) one then obtains from 
Eq.(|4]) an Efimov spectrum of tetramers: 
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The whole issue is thus to determine the values of the 
exponents s. In particular, the critical mass ratio ac(3; 1) 
corresponds to one of the exponents being equal to zero, 
the other ones remaining positive. To this end, we cal- 
culate the zero energy four-body wavefunction with no 
specific boundary condition on F{R). Then, from Eq.Q 
with i? = 0, it appears that F{R) oc R^'. The calcula- 
tion is done in momentum space, with the ansatz for the 
Fourier transform of the four-body wavefunction: 



?/'(ki,k2,k3,k4) = 



X [i?(k2,k3)+D(k3,ki)+i?(ki,k2)], (7) 

where the fermionic symmetry imposes £)(k2,ki) = 
— Z?(ki, k2), and the denominator originates from the ac- 
tion of H in Eq.([T]) written in momentum space. When H 
acts on one of the three l/|r4 — r^j singularities in Eq.(I2|), 
this produces in the right hand side of Schrodinger's 
equation a Dirac distribution 6{ri — r,) multiplied by a 
translationally invariant function of the three fermionic 
positions, which after Fourier transform gives each of the 
D[(kj)j^i^4] terms in Eq. ([?])■ Taking the Fourier trans- 
form of Eq.(|31) with F{R) cx and using a power- 
counting argument, one finds the scaling law 



D{Xki, Ak2) = A-(±^+^/2)i?(ki, k2). 
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Implementing in momentum space the contact condi- 
tions, that is the fact that 0(|ri — r4|) vanishes for = 
in Eq.([2]), gives rise to an integral equation: 
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We now use rotational invariance to impose the value 
Z S N of the total angular momentum of the four-body 
state and to restrict to a zero angular momentum along 
the quantization axis z. Then, according to Eq.©, the 
effective two-body function £'(ki,k2) has the same an- 
gular momentum /. This allows to express D in terms 
of 2Z -|- 1 unknown functions fm^ of three real variables 
only, the moduli ki and fc2 and the angle 6 € [0,7r] be- 
tween kl and k2 , with the fermionic symmetry imposing 

fi\{k2M,e) = (fci,fc2,^) Ell: 



i^(ki,k2)= J2 [Yrii,s)re^"^^'/'f,^:,iik,,b2,e). 

mi= — l 

(10) 

Here YJ™' (7, S) are the usual spherical harmonics, 7 and 
6 are the polar and azimuthal angles of the unit vec- 
tor vector along z in the direct orthonormal basis 
(ei,e2±,ei2), with ei = ki/fci,e2 = k2/fc2, e2± = 
(e2 — eicos0)/sin0 and ei2 = ei A 62/ sinO E2l- The 
action of parity k^ — > — k^ on this general ansatz is to 
multiply each term of index m; in Ea. (|10p by a factor 
(—1)™' , which allows to decouple the even mi terms (even 
parity) from the odd mi terms (odd parity). A relevant 
example, as we shall see, is the even parity channel with 
I ~ 1, where the ansatz reduces to a single term, which 
is obviously the component along z of a vectorial spinor: 



where we recall that a = M/m. Eq.® can also be ob- 
tained as the zero range limit of finite range models . 
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The last step is to use the scaling invariance of D, see 
Eq.®, setting 

/«(fcl, fc2,^) = (k! + fc2')-(^+^/2)/2(cOshx)3/2$« (X,«) 

where u = cos 9. The introduction of the logarithmic 
change of variable x = ln(fc2/fci) is motivated by Efimov 
physics, and the factor involving the hyperbolic cosine 
ensures that the final integral equation involves a Hermi- 
tian operator. The fermionic symmetry imposes 



ci>W(_:,,^) = (_l)'+i$«J:,,„) 
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which allows to restrict the unknown functions $m; to 
2: > 0. Restricting to s = iS, > 0, we finally obtain 
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The symmetrized kernel K-^J^^^ ^_^,{x,u; x' ,u') = 
Ee,e'=±i(ee')''^^^^'£™,,£',„'(e2;,w;e'a;',u') is expressed 
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FIG. 1: Minimal eigenvalues of the Hermitian operator Ms=o 
in each sector of fixed parity and angular momentum < 
I < 6, as functions of the mass ratio a = M/m. Only the 
curve for the even sector of Z = 1 crosses zero for a < 13.607, 
corresponding to the occurrence of a four-body Efimov effect 
in that sector. The other curves all remain above zero. They 
strongly overlap and are barely distinguishable at the scale of 
the figure. The dotted line is the analytical prediction A(fc — 
0, a) for the lower border of the continuum in the spectrum 
of Ms=o- The inset is a magnification. In the numerics, x and 
u were discretized with a step dx = du = 1/10, and x was 
truncated to Xmax = 20. 



in terms of the non-symmetrized one given by: 
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Here the notation V in the denominator stands for V ~ 
uu' -\- cos^Vl — u2\/l — A = e"^. A' = , is the 
angular momentum operator along x, \l,mi) is of spin I 
and angular momentum mih along z, and (f) stands for the 
azimuthal angle of the vector ka of Eq. ([9]) in the sphcrieal 
coordinates related to the basis (e2_L, ei2, ei) jlSj . 

We first look for the critical mass ratio for the 3-1-1 
fermionic problem ac(3; 1), which is the minimal value 
of a such that the integral equation Eq.(fT4|) is satisfied 
for 5 = 0. Rewriting Eq.dll]) as = Ms[$], where Ms 
is a Hermitian operator, we calculated numerically the 
minimal eigenvalues of M^^q functions of the mass ra- 
tio a, within each subspace of fixed parity and angular 
momentum <l <Q. As shown in FiglU such a mini- 
mal eigenvalue vanishes for a < 13.607 only in the even 
sector of angular momentum I ~ 1. We also unfruitfuUy 
explored I = 7,8,9, 10. We thus find that the four-body 
Efimov effect takes place only in the even sector of / = 1, 
and sets in above a critical mass ratio [l3| 

ac(3;l) ~ 13.384, (16) 

quite close to the 2 + 1 critical value ac{2; 1) ~ 13.607. 



To gain some insight on this result, we have stud- 
ied analytically an important feature of the spectrum 
of Ms=a, the lower border of its continuum. When 
x,x' — > -|-oo, which corresponds physically to having 
^2 S> ki in the function D(ki, k2), both the symmetrized 
and non-symmetrized kernels reduce to the asymptotic 
form 
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Since V is independent of x and x', this is invariant by 
translation over the x coordinates, leading to a contin- 
uous spectrum of asymptotic plane wave eigenfunctions. 
In the even sector of angular momentum I — 1, we found 
that ^'^\x, u) ~ e'^^^Vl — gives rise to an eigenfunc- 
tion in the continuous spectrum of Mis with the real 
eigenvalue A{k — S, a) |15| where 



A(fc,a) = cos 2/3 + 



(l-ik) sm[2 (3 {1 + ik)] - c.c. 
2(1-Ffc2)sin2 2/3sin(ifc7r/2) 
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In Ea. ((T5)) we have set for convenience sin 2/3 — a/{l + a) 
with /3 S [0, 7r/2[. For real k, this function A{k,a) has a 
global minimum in A: = 0. We expect that A(fc = 0, a) 
is the lower border of the continuous spectrum of Ms=o- 
Since A(0,a) exactly vanishes for the three-body critical 
mass ratio ac(2; 1) ~ 13.607, our asymptotic analysis 
amounts to uncovering the three-body problem as a limit 
^2/^1 — ^ +00 of the four-body problem. 

We tested this prediction against the numerics, plot- 
ting in FigUthe quantity A(fc = 0, a) as a function of a 
in dotted line. Except for the even sector oi I = 1, the 
minimal numerical eigenvalues are close to A(fc = 0,a); 
the fact that they are slightly above is due to a finite Xmax 
truncation effect, that indeed decreases for increasing 
2;max (not shown). This implies that the eigenfunctions 
corresponding to these minimal eigenvalues are extended, 
that is not square integrable. The numerics agrees with 
this analysis. In the even sector of / = 1, the minimal 
numerical eigenvalue is clearly below A(0,Q!), for all val- 
ues of a in Fig[TJ This indicates that the correspond- 
ing eigenvector must be a bound state of Ms^o-, with a 
square integrable cigenfunction ^'^\x,u). This is con- 
firmed by the numerics, which shows that at large x, 
$Q^^(a;, u) cx \/l — u'^e^^^ [l^. The analytical reasoning 
even predicts the link between the minimal eigenvalue 
Amin of Ms=o and the decay constant k: The plane wave 
e^''^ is analytically continuated into a decreasing expo- 
nential if one sets k = in, so that 

Amin = A(lK,a). 



(19) 



Numerically, we have successfully tested this relation for 
various values of a, and we also found that AIs=o has no 
other bound state in the even sector of Z = 1. 
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FIG. 2: In the Efimovian channel I — 1 with even parity, 
modulus of the purely imaginary Efimov exponent s — iS as 
a function of the mass ratio a = M/m. In the numerics, a;max 
ranges from 40 to 120, dx = 1/10, d9 = 7r/20. The dashed line 
results from a linear fit of |sp as a function of a in a vicinity 
of the critical value ac(3; 1), |s||t — 2.23 x (a — Oc). The 
vertical dotted line indicates the 2 + 1 critical value Qc(2; 1). 



Finally, we completed our study of the four-body Efi- 
mov effect by calculating, as a function of the mass ratio 
a, the exponent s ^ iS in the even sector of / = 1, the 
real quantity S being such that the operator Alis has a 
zero eigenvalue. The result is shovirn in FigI2] Close to 
the 2-1-1 critical mass ratio ac{2; 1) ~ 13.607, the values 
of |s| are not far from the three-boson Efimov exponent 
I So I — 1 proved to have observable effects Close to the 
3-1-1 critical mass ratio ac(3; 1), |s| varies as expected as 
{a — acY^^ (see dashed line). Low values of |s| may lead 
to extremely low Efimov tetramer binding energies: For 



an interaction of finite range 6, setting i?/ « 5 and n = 1 
in Eq.([5]), we estimate the ground state Efimov tetramer 
energy for \s\ < 1 as E^fl^ « -e-^''/\''\h^ /{2fhb'^) 
For \s\ = 0.5, taking the mass of '^He for m and a few nm 
for b gives E^^/ks in the nK range, accessible to cold 
atoms. Moreover, for a large but finite scattering length 
a, successive Efimov tetramers come in for values of a in 
geometric progression of ratio e'^/'^', so that too low val- 
ues of |s| require unrealistically large values of the scat- 
tering length. Another experimental issue is the narrow- 
ness of the mass interval. Several pairs of atomic species 
have a mass ratio in the desired interval, e.g. '^He* and 
■^^Ca {a ~ 13.58), and with exotic species, ^"'^B and ^'^^Sm 
(a ~ 13.53), '^Li and ^^Mo {a ~ 13.53). A more flexi- 
ble solution is to start with usual atomic species having a 
slightly off mass ratio, such as •^Hc* and {a ~ 13.25), 
and to use a weak optical lattice to finely tune the effec- 
tive mass of one of the species [l3| ■ 

To conclude, in the zero range model at unitarity, we 
studied the interaction of three same spin state fermions 
of mass Af with another particle of mass m. For M/m < 
13.384, no Efimov effect was found. Over the interval 
13.384 < M/m < 13.607, remarkably a purely four-body 
Efimov effect takes place, in the sector of even parity and 
angular momentum ^ = 1, that may be observed with a 
dedicated cold atom experiment. For M/m > 13.607, 
the three-body Efimov effect sets in, and the zero range 
model has to be supplemented by three-body contact con- 
ditions that break its separability. The intriguing ques- 
tion of wether the Efimov tetramers then survive as res- 
onances, decaying in a trimer plus a free atom, is left for 
the future. F. Werner is warmly thanked for discussions. 
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